The "unusual " isotope shift in high-temperature superconductors 
can be explained by the usual theory of the electron-phonon interaction 
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We show that recent ARPES resuhs on the "unusual" oxygen isotope shift in the real part of the 
self-energy in the optimally doped _Bi2212 samples Q can be qualitatively (and semi-quantitatively) 
explained by the theory of the electron-phonon interaction (EPI) elaborated few decades ago. How- 
ever, for a quantitative analysis of the ARPES spectra it is necessary to know the momentum 
dependence of the EPI, the Coulomb contribution at high energies and the background due to 
impurities and defects. 

PACS numbers: 



Recently a very interesting paper was reported in the 
"Nature" with the title "An unusual isotope shift in 
high-temperature superconductors" by Gwon et al 0|, 
where the angle-resolved photoemission (ARPES) spec- 
tra have been investigated in the optimally doped Bi2212 
samples for different stages of the oxygen isotope substi- 



O. It was shown that all energy 



tution O 
distribution curves (EDS) demonstrate a rather small but 
nonzero isotope effect. Each of these curves shows a peak, 
the position of this peak is affected by the above men- 
tioned isotope substitution. All isotope effects have max- 
imal values for frequencies of the order of 100-300 meV 
and vanish for larger energies. 

In this short comment we demonstrate that the main 
part of ARPES results in Ref.Jj] can be easily un- 
derstood, and even semi-quantitatively explained, in 
the framework of the standard Migdal-Eliashberg model 
|M 1^ Q for the electron-phonon interaction (EPI). For 
that purpose we discuss below the quasiparticle self- 
energy S(w) due to the EPI in the simple model in 
which the quasiparticles interact with a dispersionless 
optical phonon (Einstein model) , which has been studied 
at length in Ref. ^4] niany years ago. In this analitically 
solvable model the Eliashberg spectral function has the 
form a'^F{uj) = Xujq6{uj — ojq) Below we shall discuss the 
energy and isotope dependence of ReT,{Lu) for the quasi- 
particle momenta in the nodal direction (0,0) — (tt, tt), 
since in this case there are no additional effects on S(ijj) 
due to the superconducting gap and the pseudogap. The 
obtained results in the simple EPI model will be com- 
pared with those of the ARPES spectra in Q . 

The real part of S(w) in the simple EPI model has the 
following form Q 



i?eE(c.) = -^ln 
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(1) 



where A = g'^N{0/ujQ is the EPI coupling constant, g is 
the matrix element of the EPI, N{Q) is a quasiparticle 
density of states on the Fermi level and loq is the optical 
phonon energy. We would like to mention here that the 



value ReTtluj) is negative, while | i?eS](a;) |was plotted in 
the inset of Fig. 1 in Ref. 1]. 

From Eq.(^ it is seen that |i?eS(Li;)| in the simple 
model is logarithmically singular at w ~ wq- It is well 
known that this singularity is smoothed by (i) the real- 
istic phonon spectrum and the Eliashberg spectral func- 
tion, and (ii) by the nonzero temperature effects T ^ 

. This will be studied elsewhere for spectral dunctions 
related to high-temperature superconductors. 

The position of the peak in |i?eS(a;)| depends on iso- 
tope substitutions because the phonon energies depend 
on atomic masses 



(2) 



Here 7 is a force constant and M is an effective mass of a 
phonon mode. The high-energy optical phonon modes in 
high- Tc superconductors are mainly related to the oxygen 
motion. This means that the value M in Eq.(|21) is an 
oxygen mass. Ea. (|ll2|l show the red shift of the peak 
position for the heavier isotope. The value of this shift 
is equal 

Auj = cjo(^^O) - ujoi^^O) ~ QMljoC^O) (3) 

whose order of magnitude is in a good agreement with 
the observation in Ref. . 

Now we consider |i?eE(w)| at low energies (oj <C ujq) 
where Eq.(^ gives a linear behavior 



\ReY.{Lu) \ ~ Xuj. 



(4) 



It is important to stress that in metals A does not depend 
on the mass of vibrating atoms in the adiabatic and har- 
monic approximation, as it was first shown in j^]. We 
stress that this result was observed in ARPES spectra of 
Ref. 1], where the slope of |i?eI](Li;)| sd lu <^ ujq is isotope 
independent. It means, that the high-temperature super- 
conductors are not in the so called nonadiabatic regime 
of the polaron formation. The expression Eq. JQ) can be 
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expanded up to the next power of uj /ljq where 

|i?eS(w)| ^ Aw(l - — ). (5) 

The smah deviation (from the hnear dependence) of 
i?eS(a;)l is negative and is larger for the heavier iso- 
tope. This result is also in a good agreement with the 
ARPES results shown in the inset of Fig. (lb) in Ref.0. 

Let us now consider |i?el]((x')| at high energies uj ^ wq. 
In this case one has 

|i?eE(w)| cc Aa;o( — ). (6) 

UJ 

This expression shows that |i?eS](w)| l/cj, i.e. it falls 
off slowly by increasing uj. Moreover |i?eE(w)| ^ 
i.e. it depends on the isotope mass. It means that the 
value of the isotope effect is more pronounced at high 
energies then at low. This result is also in a qualita- 
tive agreement with the observed ARPES data in Ref . T| . 
Note, the same analitical behavior should hold for more 
realistic spectral functions, than for the Einstein model), 
for w ^ LOmax, where uJmax is the maximal phonon en- 
ergy. 

The simple theory presented here cannot explain the 
absolute value of the isotope effect at very large w, which 
was observed in Ref. , since we do not know the contri- 
bution of the Coulomb interaction to |i?el](w)| at these 
energies. However, we can use the difference between 
the experimental values of |i?eS([x))| for ^^O and ^^O 

at rather high energy uj — 0.2eV, i.e. 6 \ReT,{uj)\ieo = 

16q 18q 

\ReT,{uj)\ — |i?el]((jj)| , in order to estimate the EPI 
coupling constant A(;j) (in the nodal direction) in the 
above simple model. In absence of much more detailed 
experimental details in this energy region, we estimate 



the experimental value of 5 \ReT,{w) \ \]_Iq from the small 
inset in Fig.l of Ref. Q to be |i?eS(cj)| |ieg ~ (6 - 10) 
meV, which combined with Eq.lQ, gives A(;i) ~ (2 — 3). 
This value is at least by factor two larger than the cou- 
pling A(;) ^ 1 extracted from the slope of |i?eE(w)| at 
low w < tJo in the ARPES data Q , [3 • Although this 
discrepancy in coupling constants A(;) andA(;i) is not too 
large, it requires more accurate experiments and theoret- 
ical analysis in the high-energy region. 

We stress that for a quantitative analysis of the 
ARPES data in high-temperature superconductors it is 
necessary to know (1) the EPI spectral function and its 
momentum dependence, (2) the contribution from the 
Coulomb interaction, (3) the background scattering, etc. 
For instance, the point (1) is important in order to ex- 
plain the so-called non-shift ARPES puzzle, where the 
quasiparticle kink in the nodal direction at 70 mcV is un- 
shifted in the superconducting state, while the one near 
the anti-nodal point (at 40 meV) is shifted. This non- 
shift puzzle was explained in Q by the existence of the 
forward scattering peak in the EPI and impurity scatter- 
ing. 

In conclusion, the theory of the electron-phonon in- 
teraction is capable in explaining the recent ""unusual" 
isotope shift of the real part of the self-energy of high- 
temperature superconductors 1], without invoking exotic 
mechanisms for the quasiparticle interaction. However, 
in spite the fact that recent ARPES data 0, 0, H, 
10] favor the EPI as the pairing mechanism in the high- 
temperature superconductors for a quantitative analy- 
sis the knowledge of much more microscopic details are 
needed. 
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